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Preface 


The  purpose  of  this  study  was  to  measure  the  total 
hemispherical  emittance  of  the  fuel  cladding  for  the  SP-100 
space  nuclear  reactor  core.  The  reason  for  measuring  this 
property  was  to  provide  a  more  accurate  emittance  value  for 
computer  modeling  of  the  consequences  of  a  loss-of-coolant 
accident  occurring  in  space.  Although  the  emittance  of  Nb- 
1%  Zr  had  been  measured  by  others,  there  remained  some 
degree  of  uncertainty  concerning  the  possible  effects  of 
cladding  geometry  and  core  environment  on  the  alloy's  ther¬ 
mal  radiation  properties.  This  study  was  designed  to  pro¬ 
vide  a  basis  for  conducting  emittance  research  on  the  clad¬ 
ding  by  attempting  to  measure  the  emittance  of  an  actual 
sample  in  a  vacuum. 

The  method  and  apparatus  for  calculating  emittance  were 
based  on  a  standard  test  method  (ASTM  C  835-82,  Reapproved 
1988) .  Although  emittance  was  reported  over  the  temperature 
range  of  about  900  to  1300  K,  the  more  interesting  region 
toward  the  upper  core  operating  limit  of  2000  K  was  not 
explored  due  to  equipment  limitations.  In  addition,  the 
specimen  was  affected  by  gradual  oxidation  resulting  from 
heating  under  a  less-than-adequate  vacuum,  a  condition  which 
hampered  statistical  error  analysis.  Despite  such  limita¬ 
tions,  the  method  provided  conclusive  results.  With  some 
procedural  modifications  and  improved  equipment,  it  appears 
promising  for  continued  research. 
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Abstract 


The  purpose  of  this  study  was  to  measure  the  total 
hemispherical  emittance  of  Nb-1%  Zr  fuel  cladding  for  the 
SP-100  space  nuclear  reactor  in  the  operational  temperature 
range.  An  apparatus  was  constructed  based  primarily  on  ASTM 
C  835-82,  Standard  Test  Method  for  Total  Hemispherical  Emit¬ 
tance  of  Surfaces  from  20  to  1400*  C  (Reapproved  1988) .  The 
test  method  required  heating  the  specimen  in  a  vacuum  until 
achieving  thermodynamic  equilibrium  between  the  specimen  and 
its  surroundings.  This  was  accomplished  by  passing  an 
electric  current  through  the  specimen  and  then  allowing  its 
surface  temperature  to  stabilize.  Measurements  were  made  of 
the  surface  temperature  and  the  electric  power  dissipated 
within  a  predetermined  region.  The  emittance  was  then 
calculated  using  the  Stefan-Boltzmann  Law  and  the  assumption 
that  the  power  was  dissipated  solely  by  thermal  radiation. 

An  oxidized  stainless  steel  specimen  of  known  emittance  was 
used  for  verification. 

For  the  virgin  (unoxidized)  specimen,  measured  emit¬ 
tance  ranged  from  0.159  ±  1.34%  at  913  K  to  0.200  ±  1.25%  at 
1091  K.  The  highest  emittance  value  obtained  was  0.339  ± 
1.63%  at  1269  K  for  the  specimen  after  visible  surface 
oxidation  was  noted.  The  available  vacuum  was  found  to  be 
inadequate  for  preventing  oxidation  of  the  heated  Nb-1%  Zr, 
causing  a  steady  increase  in  the  emittance  value  of  the 
specimen  from  one  experiment  run  to  the  next. 
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TOTAL  HEMISPHERICAL  EMITTANCE  OF  NIOBIUM-1%  ZIRCONIUM 
FUEL  CLADDING  FOR  THE  SP-100  SPACE  REACTOR 


X-  Introduction 


Background 

The  SP-lOO  Program,  as  described  by  J.  Angelo  and  D. 
Buden  in  Space  Nuclear  Power,  is  a  venture  begun  in  1983  by 
the  Department  of  Defense,  the  Department  of  Energy,  and  the 
National  Aeronautics  and  Space  Administration  to  "jointly 
develop  space  reactor  technology”  (3:232).  Goals  for  the 
program  include  a  net  electric  power  output  of  100  kW,  a 
mass  constraint  of  3,000  kg,  a  maximum  length  of  6.1  m,  and 
safe  operation  (3:233).  Throughout  the  entire  project,  the 
design  guidelines  have  included  a  requirement  for  ensuring 
full  consideration  of  aerospace  nuclear  safety.  This  study 
concerns  one  safety  aspect:  the  question  of  how  well  the  Nb- 
1%  Zr  cladding  can  radiate  decay  heat  from  the  core  follow¬ 
ing  the  loss  of  its  primary  coolant. 

The  coolant,  liquid  lithium,  circulates  under  pressure 
around  the  fuel  pins  (uranium  nitride  pellets  clad  in  tubu¬ 
lar  sleeves  of  Nb-1%  Zr)  (3:234).  Normally,  the  coolant 
leaves  the  reactor  core  at  temperatures  between  1131  K  and 
1375  K  (8:1).  Since  a  space  environment  subjects  the  reac¬ 
tor  to  bombardment  by  meteors  and  orbiting  space  debris,  the 
cooling  system  is  vulnerable  to  puncture.  Such  an  event 
would  evacuate  the  space  around  the  pins  and  is  known  as  a 
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"loss-of“COolant  accident”  or  LOCA.  Without  the  presence  of 
the  coolant,  the  danger  exists  that  the  maximuin  core  temper¬ 
ature  of  2000  K  could  be  exceeded  (4:575). 

Following  a  LOCA  and  an  emergency  reactor  shut  down, 
the  core  must  survive  the  heat  accumulating  from  decaying 
fission  products.  In  order  to  control  the  temperature  of 
the  core  and  prevent  structural  failure,  the  decay  heat, 
initially  produced  at  the  rate  of  about  seven  percent  full 
power,  would  need  to  be  removed.  One  solution  uses  an 
auxiliary  supply  of  lithium  that  is  isolated  from  the  prima¬ 
ry  reservoir  and  connected  to  tubes  interspersed  among  the 
fuel  pins.  For  maximum  effectiveness,  these  tubes,  called 
cold  fingers  or  bayonets,  would  be  arranged  within  the  core 
so  that  each  can  dispose  of  heat  from  several  nearby  fuel 
pins.  In  addition,  no  pin  must  be  allowed  to  exceed  the 
2000  K  temperature  limit. 

Robbins  (14)  used  T.  Kao's  Thermal  System  Analysis 
Program  (TSAP)  code  (10),  in  modeling  the  proper  placement 
of  the  bayonets,  based  on  assumed  values  of  emittance.  The 
important  trend  demonstrated  by  TSAP  is  that  higher  emit¬ 
tance  values  for  the  fuel  cladding  corresponded  to  a  rela¬ 
tively  larger  zone  of  influence  for  any  one  bayonet.  It  was 
thus  shown  that  the  total  number  of  required  bayonets  de¬ 
creased  as  the  emittance  increased. 
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Definition  q£.  Einittance 

Emittance,  often  called  thermal  emittance,  is  the 
property  of  a  material  that  relates  the  amount  of  energy 
radiated  to  that  of  a  black  body  emitter  at  the  same  temper- 
ature.  A  point  of  possible  confusion  exists  in  that  emit¬ 
tance  is  often  used  interchangeably  with  "emissivity" . 
However,  one  system  of  nomenclature  in  use  reserves  the 
"ivity"  suffix  to  intensive  properties  describing  a  material 
in  general.  Extensive  properties,  such  as  emittance,  would 
be  used  to  describe  a  specimen  of  the  material  (7:3-4).  For 
instance,  emissivity  would  apply  to  the  alloy  Nb-1%  Zr, 
while  emittance  might  be  used  in  describing  the  surface  of  a 
12-inch-long,  thin-walled,  hollow  cylinder  of  polished  Nb-1% 
Zr. 

According  to  the  Standard  Test  Method  for  Total  Hemi¬ 
spherical  Emittance  of  Surfaces  from  20  to  1400  "C,  set 
forth  by  the  American  Society  for  Testing  and  Materials 
(ASTM  C  835-82,  Reapproved  1988),  emittance  can  be  charac¬ 
terized  into  normal  vs  hemispherical  and  spectral  vs  total . 
Normal  emittance  corresponds  to  energy  radiating  normal  from 
the  surface,  while  hemispherical  emittance  relates  to  energy 
radiating  at  all  angles  (2:1).  Some  applications  require 
characterization  of  the  energy  spectrum.  In  those  cases, 
spectral  emittance  would  be  used  to  describe  an  object's 
radiation  properties  as  they  apply  to  a  specific  spectral 
band.  Accordingly,  total  emittance  would  include  all  wave¬ 
lengths.  Total  hemispherical  emittance  is  used  here,  be- 
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cause  the  interest  lies  in  all  energy  being  radiated  at  all 
angles  from  the  bare  fuel  pin  cladding.  Figure  1  shows 
three  views  of  radiation  from  a  point  on  the  cylinder  sur¬ 
face  and  demonstrates  the  hemispherical  nature  of  such 
radiation. 


Figure  1.  Hemispherical  Radiation  from  a  Point  on  the 

Surface  of  a  Cylinder 


Problem 

This  study  is  an  attempt  to  adapt  the  method  prescribed 
for  thin  strips  by  ASTM  C  835-82  for  use  on  a  thin  cylin¬ 
drical  test  specimen.  A  working  apparatus  conforming  to  the 
standard  was  located  at  ORNL,  and  a  visit  was  made  to  in¬ 
spect  this  apparatus  and  the  experimental  procedure.  Since 
a  long,  thin-walled  cylinder  could  be  considered  a  thin 
metal  strip  bent  to  form  a  continuous  surface,  it  was  deter¬ 
mined  that  the  ORNL  design  would  work.  Several  assumptions 
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were  required  in  order  to  adapt  the  design  to  this  study 
while  maintaining  compliance  with  the  ASTM  standard;  the 
assumptions  are  discussed  separately  in  section  I.  Figure  2 
shows  an  overview  of  the  apparatus. 


Figure  2.  Overview  of  Apparatus 


This  study  provides  values  for  the  total  hemispherical 
emittance  of  a  sample  specimen  of  SP-100  specification 
cladding  obtained  from  Los  Alamos  National  Laboratory 
(LANL) .  The  sample  was  approximately  30  cm  long  with  an 
outside  diameter  of  0.76035  cm.  Type  S  thermocouples  were 
first  attached  to  the  sample.  Also  serving  as  supports, 
copper  clamps  were  used  inside  the  bell  jar  to  connect  power 
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lines  from  the  DC  source  to  the  sample.  The  clamps  and 
cables  were  required  to  be  large  enough  to  handle  currents 
up  to  180  amps  for  heating  the  sample. 

In  order  to  obtain  a  value  for  emittance,  the  voltage 
drop  across  a  predetermined  section  of  the  sample  was  mea¬ 
sured,  along  with  the  current  flowing  through  the  sample  and 
its  surface  temperature.  A  discussion  appears  in  Section  II 
concerning  how  these  physical  values  were  used  in  calculat¬ 
ing  emittance. 

Scope 

Senor  and  others  correlated  existing  data  sets  and 
demonstrated  a  temperature  dependence  for  total  hemispheri¬ 
cal  emittance  of  the  subject  alloy  (15:206).  Ideally,  this 
dependence  should  be  characterized  for  cladding  surface 
temperatures  up  to  2000  K,  the  upper  operating  limit  for  the 
SP-100  core.  However,  this  study  limited  the  temperature  to 
1268  K  maximum  due  to  power  supply  constraints.  The  appara¬ 
tus  initially  included  a  75  amp  DC  power  supply,  replaced  by 
a  180  A  supply  for  the  final  series  of  measurements. 

Another  limiting  factor  was  the  availability  of  suit¬ 
able  vacuum  equipment.  To  avoid  convective  loss,  the  vacuum 
chamber  required  a  pressure  of  1.3  mPa  or  less  (2:3). 
Otherwise,  convective  loss  would  drive  down  the  surface 
temperature  and  increase  the  apparent  emittance.  In  addi¬ 
tion,  the  vacuum  prescribed  by  the  ASTM  standard  was  insuf¬ 
ficient  for  Nb-1%  Zr,  which  exhibited  surface  oxidation 
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problems  at  pressures  as  low  as  0.53  mPa.  The  steady  oxida¬ 
tion  of  the  sample  gradually  increased  the  emittance  over 
time. 

Assumptions 

The  derivation  used  in  this  method  to  calculate  emit¬ 
tance  required  the  consideration  of  several  basic  assump¬ 
tions.  The  first  three  were  based  on  the  ASTM  standard. 

They  are: 

1.  The  enclosure  is  a  blackbody  emitter/absorber  at  a 
uniform  temperature.  The  ASTM  standard  permits  the  use  of 
an  uncoated  glass  bell  jar  because  "it  absorbs  or  transmits 
most  of  the  incident  infrared  radiation"  (2:9). 

2.  The  sample's  total  hemispherical  emittance  equals  its 
absorptance.  This  assumption  refers  to  the  exchange  of 
diffuse  blackbody  radiation  between  the  sample  and  its 
surroundings . 

3.  There  is  no  heat  loss  from  the  test  section  by  convec¬ 
tion  or  conduction  (2:2). 

Section  II  contains  a  discussion  of  these  assumptions 
in  greater  detail,  as  well  as  a  fourth  assumption  required 
due  to  the  geometry  involved. 

4.  The  dimensions  of  the  cylindrical  specimen  are  suffi¬ 
cient  to  prevent  more  than  a  negligible  net  gain  or  loss  of 
radiative  energy  at  the  inner  surface  of  the  test  section. 
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Apprgftgh 


To  establish  a  basis  for  this  study.  Section  II  will 
begin  with  the  development  of  the  theory  involved  in  emit- 
tance  measurement.  A  brief  introduction  explains  the  rele¬ 
vance  of  emittance  to  radiative  heat  transfer.  Following 
the  introduction,  all  relevant  heat  transfer  mechanisms  are 
applied  to  the  problem. 

Following  the  theory.  Section  III  covers  the  design  and 
construction  of  the  apparatus.  The  assumptions  made  earlier 
will  be  treated  in  greater  detail  to  explain  how  they 
affected  the  design.  Section  IV  explains  the  experimental 
procedure.  Included  in  this  section  are  the  procedures  for 
operating  the  apparatus  and  for  gathering  data.  The  effect 
of  design  constraints  on  data  collection  will  be  discussed. 

Section  V  contains  the  results  of  the  study,  along  with 
an  error  analysis  and  the  conclusions  reached.  Included  in 
this  section  is  a  comparison  of  the  AFIT  results  with  those 
from  ORNL  for  the  ORNL  oxidized  stainless  steel  emittance 
standard.  The  Nb-1%  Zr  results  are  compared  with  Senor's 
correlation  mentioned  earlier.  Finally,  Section  VI  summa¬ 
rizes  the  conclusions  made  and  contains  recommendations  for 
further  research  on  Nb-1%  Zr  fuel  cladding  emittance. 
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Introduction 

This  section  contains  the  theory  underlying  the  assump¬ 
tions  made  when  measuring  total  hemispherical  emlttance. 
First,  the  Stefan-Boltzmann  Law  will  be  used  to  highlight 
the  difference  between  a  blackbody  radiator  and  less-than- 
Ideal  radiators,  reemphasizing  the  relevance  of  thermal 
radiation  to  the  study.  Next,  the  heat  transfer  mechanisms 
present  In  a  generic  emlttance  laboratory  experiment  will  be 
discussed  as  a  basis  for  the  method  used.  Finally,  the 
theory  will  be  adapted  to  the  particular  case  at  hand,  a 
hollow  cylindrical  radiator  of  a  refractory  alloy. 

Stefan-Boltzmann  Law 

The  Stefan-Boltzmann  Law  states  that  the  power  per  unit 
area  radiated  by  a  "blackbody"  or  Ideal  radiator  Is  propor¬ 
tional  to  the  fourth  power  of  Its  absolute  surface  tempera¬ 
ture.  This  relationship  Is  given  by 

(1) 


where 

P,  =  power  radiated  per  unit  area 
a  =  Stefan-Boltzmann  Constant 
T  ==  absolute  temperature 

To  account  for  the  loss  of  radiation  efficiency  when  the 
radiator  Is  less  than  Ideal,  one  uses  a  factor  known  as 
emlttance,  or  the  ratio  of  the  power  per  unit  area  radiated 
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by  the  non-blackbody  to  that  radiated  by  a  blackbody  at  the 
same  temperature.  Eg  (1)  thus  becomes 

=  ear*  (2) 

where  c  is  the  emlttance,  a  positive  number  less  than  one 
for  all  non-ideal  radiators. 

Now  that  the  background  has  been  established,  the 
remainder  of  Section  II  will  be  used  to  derive  the  rela¬ 
tionships  used  in  this  study. 

Heat  Flow  In  n  Solid  Cylinder  Conducting  Mgtal 

Consider  a  solid,  conducting  metallic  cylinder  with  DC 
current  applied  to  the  ends  via  electrodes.  Due  to  its  own 
resistivity,  the  cylinder  will  begin  to  produce  thermal 
energy,  and  several  heat  transfer  processes  will  attempt  to 
return  the  cylinder  and  its  surroundings  to  thermal  equilib¬ 
rium.  Besides  radiation,  these  processes  include  conduc¬ 
tion,  convection,  absorption,  and  the  storage  of  internal 
energy  (1:294). 

A  simple  expression  for  the  transport  processes  men¬ 
tioned  above  per  unit  length  is 

Pd‘  Pr-  P.*  Pai*  Pcv^  Ps  (3) 
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where 


P(,  =  electric  power  dissipated  per  unit  length 

Pp  =  power  radiated  away  per  unit  length 

Pj,  =  power  absorbed  from  surroundings  per  unit  length 

Pcd  =  power  lost  through  conduction  per  unit  length 

Pcv  =  power  lost  through  convection  per  unit  length 

Pj  =  power  lost  to  internal  energy  per  unit  length 

These  terms  shall  be  treated  one  at  a  time  in  the  general 

case  of  the  solid  cylinder,  followed  by  an  adaption  to  the 

specific  case  of  the  hollow  cylindrical  cladding. 

First,  the  rate  of  increase  in  internal  energy  stored 

in  the  cylinder  shall  be  examined.  This  heat  transfer 

mechanism  is  given  as 


Ps  =  DAc-^ 
^  dt 


(4) 


where 

Pg  =  power  stored  per  unit  length  in  W  m  ’ 

D  =  density  in  kg  m'^ 

A  =  cross  sectional  area  in  m^ 
c  =  specific  heat  in  J  kg  ’  k  ’ 

T  =  absolute  temperature  in  K 
t  =  time 

In  order  for  this  term  to  have  any  effect,  the  temperature 
must  be  changing  with  respect  to  time.  At  this  point,  a 
boundary  condition  shall  be  imposed  on  the  problem;  the 
derivative  of  temperature  with  respect  to  time  must  be  0,  or 
in  other  words,  a  steady-state  condition  must  be  present. 
Thus,  the  Pg  term  is  eliminated. 

The  next  term  to  be  discussed  shall  be  the  conduction 
loss  term.  In  order  for  heat  conduction  to  occur,  a  tem¬ 
perature  gradient  must  exist.  Therefore,  if  it  can  be  shown 
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that  no  temperature  gradient  exists  along  a  given  coordi¬ 
nate,  it  can  be  said  that  the  conduction  loss  term  contains 
no  component  associated  with  that  coordinate.  To  assume 
that  no  net  conduction  loss  occurs  from  a  particular  zone  of 
interest  in  a  cylindrical  coordinate  system,  then,  one  must 
show  that  no  temperature  gradient  exists  along  the  radial, 
azimuthal,  or  axial  coordinates.  The  following  section 
contains  a  treatment  of  this  problem  as  it  applies  to  the 
specimen  used  in  emlttance  measurement. 

Ordinarily,  one  might  use  LaPlace's  equation  to  de¬ 
scribe  the  steady-state  temperature  distribution  in  a 
cylinder.  In  cylindrical  coordinates,  it  is  given  as 


1  .  0 
p  3p  V  3p/  p*  30^  Bz^ 


(5) 


where  p,  0,  and  z  are  the  radial,  angular,  and  longitudinal 
coordinates  respectively.  The  basic  solutions  of  interest 
for  T  are 


T  = 


(icp)  sin(n0)exp(-icz) 
(icp)  cos  (j36>  exp  (-1:z) 


1 


(6) 


where  k  (real)  is  a  zero  of  the  Bessel  function  J„  and  the 
separation  constant  n  is  an  integer.  However,  for  the 
purposes  of  this  study,  the  dimensions  and  location  of  the 
zone  of  interest  will  be  chosen  such  that  the  temperature 
within  the  zone  has  no  dependence  on  p,  6,  or  z. 
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FIqw  In  n  Thin/  HqIIQW  CYiin<;i^r  Ol  conducting  Mgt»l 


Consider  a  30  cm  long  section  of  SP-100  fuel  cladding, 
a  thin-walled  hollow  cylinder  with  (1)  a  wall  thickness  very 
small  in  relation  to  the  cylinder's  diameter  and  length,  and 
(2)  a  diameter  very  small  in  relation  to  the  cylinder's 
length. 

The  discussion  will  begin  with  heat  conduction  in  the 
radial  direction.  The  analogy  to  the  ASTM  emittance  calcu¬ 
lation  standard  is  the  heat  transfer  through  the  thickness 
of  the  flat  metal  strip.  The  ASTM  standard  assumes  a  sample 
thickness  too  insignificant  to  allow  a  temperature  gradient 
to  exist  between  the  two  flat  sides,  implying  no  heat  con¬ 
duction  normal  to  the  plane  of  the  specimen.  In  other 
words,  all  heat  is  considered  to  be  generated  on  the  two 
flat  surfaces  of  the  strip.  Taking  a  similar  approach  with 
the  cylindrical  specimen,  one  can  assume  a  wall  thickness 
such  that  all  heat  is  generated  on  the  inner  and  outer 
surfaces  with  no  temperature  gradient  between  them. 

Next,  the  angular  component  of  heat  conduction  will  be 
examined.  If  one  assumes  that  the  cylinder  is  homogeneous 
and  uniform,  there  is  no  basis  for  a  temperature  gradient  to 
exist  as  a  function  of  0.  Thus,  no  angular  component  of 
heat  conduction  exists. 

Now  that  the  radial  and  angular  dependencies  of  the 
temperature  have  been  eliminated,  only  the  question  of  a 
longitudinal  dependence  remains.  At  this  point,  an  assump¬ 
tion  is  made  that  the  cylinder  is  very  long  in  relation  to 
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its  diameter.  By  measuring  the  surface  temperature  of  the 
cylinder,  one  can  find  a  region  (Figure  3)  near  the  middle 
such  that  the  axial  temperature  gradient  (caused  by  conduce 
tion  out  the  electrodes)  is  essentially  zero  (2:9).  The 
zone  of  interest,  or  test  zone,  for  this  study  shall  be 
chosen  so  that  its  boundaries  are  within  this  region  lacking 
an  axial  temperature  gradient. 


Cond. 
Loss  ^ 


Test  Region 


Power 

Source 


Relative 

Temperature 


Sample 


■  Cond. 
^  Loss 


Figure  3.  Axial  Temperature  Distribution  Along  an  Electri¬ 
cally  Heated  Metallic  Cylinder 


Since  it  has  been  shown  that  no  temperature  gradient 


exists  across  the  boundaries  of  the  test  zone,  it  can  there 
fore  be  said  that  there  can  be  no  net  heat  conduction  loss 
from  the  zone.  Thus,  the  term  in  Eq  (3)  is  eliminated. 


Before  moving  on  to  convective  loss,  it  is  appropriate 
to  mention  the  heat  conductive  loss  caused  by  attaching 
thermocouples  to  the  sample  within  the  test  zone.  Such  loss 
creates  a  temperature  depression  at  the  weld.  Since  thermo¬ 
couple  wires  of  0.003"  diameter  are  welded  to  the  specimen's 
surface  in  the  zone  of  interest,  any  conduction  loss  out 
these  wires  would  lower  the  actual  surface  temperature  in 
the  weld  region.  Using  the  ASTM  standard,  it  is  assumed 
that  any  loss  through  such  small  diameter  wire  would  be 
negligible.  In  fact,  the  ASTM  standard  states  that  "the 
heat  conducted  out  through  the  thermocouple  leads  is  usually 
neglected"  (2:9). 

Next,  the  convection  term  in  Eg  (3)  shall  be  con¬ 
sidered.  Since  the  study  attempted  to  simulate  a  space 
environment,  a  vacuum  bell  jar  was  used.  The  convection 
heat  transfer  coefficient  is  zero  in  a  vacuum,  immediately 
eliminating  the  convection  term. 

Eg  (3)  has  now  been  reduced  to 

Pd  =  Pr-Pa  (7) 

Pp,  the  power  radiated  per  unit  length,  can  be  expressed  as 

Pr^^hA:,oT*  (8) 


where 

e„  =  total  hemispherical  emittance 

Al  =  radiating  area  per  unit  length  of  specimen 


15 


A  reciprocal  heat  transfer  process  occurs,  wherein  the 
radiator  can  absorb  heat  from  its  surroundings,  assumed  here 
as  a  black  body  ed>sorber/ radiator.  This  process  is  given  by 

(9) 


where 

s  power  absorbed  per  unit  length  of  specimen 
Cif  s  total  hemispherical  edssorptance  of  the  specimen 
Tg  =  absolute  temperature  of  the  surroundings,  includ¬ 
ing  the  jar  and  the  ambient  air 

The  electric  power  dissipated  per  unit  length,  Pg,  is  simply 

the  current  flowing  through  the  specimen  times  the  voltage 

drop  per  unit  specimen  length. 

Substituting  Eqs  (8)  and  (9)  into  Eg  (7)  yields 

AjOTS  (10) 

By  multiplying  the  equation  through  by  unit  length  and 
dividing  through  by  aA,  one  obtains  the  relationship 

Pp  =  Aa{€g  T*-aj,  To)  (11) 

Eq  (11)  relates  the  radiative  heat  transfer  into  and 
out  of  a  specified  area  with  the  electric  power  dissipated 
within  that  area.  The  right  side  can  be  regrouped  to  obtain 

Ao^r*-- ^2^1  (12) 
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According  to  Abbott,  the  term  in  Eq  (12)  containing  Tq 
is  often  ignored,  accounting  for  an  error  of  less  than  1% 
when  T  is  at  least  3  times  T,,  (1:295).  Attempted  here 
initially,  the  approximation  introduced  errors  of  up  to  3% 
in  this  study,  and  was  abandoned.  The  ASTM  standard  allows 
the  use  of  an  uncoated  glass  bell  jar  by  assuming  that  the 
energy  incident  on  the  bell  jar  is  mostly  either  edssorbed  or 
transmitted.  Therefore,  Tg  was  assumed  to  be  the  ambient 
air  temperature  of  295  K.  In  addition,  the  ratio  of  edssorp- 
tion  to  emittance  was  deemed  unity  by  assximption  2  on  page 
7.  Solving  for  then,  gives 
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H 


Aa{T*  -  T^) 


(13) 


the  expression  used  to  calculate  total  hemispherical  emit¬ 
tance  in  this  study. 

One  final  assumption  was  made  concerning  radiation.  As 
stated  above,  the  cylinder  was  assumed  to  have  a  small 
diameter  relative  to  its  length.  The  purpose  for  this  was 
to  eliminate  net  radiation  gain  or  loss  at  the  zone  of 
interest's  inside  surface.  It  is  therefore  assumed  that  a 
narrow,  long,  hollow  cylinder  would  facilitate  a  condition 
of  local  thermodynamic  equilibrium  inside  the  tube  away  from 
the  ends.  See  the  recommendations  section  for  comments 
concerning  this  aspect  of  the  problem. 
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III.  Apparatus 


Ths.  fia&iss 

As  shown  in  Figure  1,  the  apparatus  for  measuring  total 
hemispherical  emittance  was  comprised  of  a  vacuum  system, 
bell  jar,  DC  power  supply,  specimen  holder,  precision  volt¬ 
meter,  thermocouples,  and  temperature  reference  bath.  A 
multiposition  switch  box  allowed  data  to  be  taken  without 
disturbing  any  leads  or  connections. 

Vacuum  System 

A  six-inch  water-cooled  diffusion  pump  was  used  in  the 
vacuxim  system,  backed  by  a  mechanical  roughing  pump.  To 
begin  with,  two  problems  associated  with  this  system  had  to 
be  solved.  First,  to  prevent  contamination  of  the  bell  jar 
and  contents  from  any  backstreaming  pump  oil,  a  cold  trap 
was  used.  It  was  maintained  at  around  -8*  C  using  a  stan¬ 
dard  refrigerator  compressor  and  R-12  coolant.  The  cold 
trap  added  its  own  problem:  an  insufficient  load  to  the 
compressor  to  prevent  freezing  of  the  coolant  return  line. 
Electrical  heat  tape  was  applied  to  the  line  to  prevent 
freezing. 

An  ionization  gauge  was  used  for  monitoring  the  vacuum 
in  the  bell  jar.  It  could  be  isolated  from  the  pumping 
system,  as  required,  to  perform  static  pressure  checks. 
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Bell  Jar 

A  quartz  bell  jar  with  a  wall  thickness  of  one-eighth 
inch  was  used  to  house  the  specimen.  For  safety  reasons,  a 
protective  metal  screen  was  installed  on  the  outside,  sur¬ 
rounding  Che  jar.  The  inner  walls  were  not  coated,  and  no 
water  cooling  was  employed. 

The  requirement  to  keep  the  recommended  100  K  differ¬ 
ence  between  the  bell  jar  and  the  specimen  was  accomplished 
at  ORNL  by  using  a  water-cooled,  coated  metallic  bell  jar 
(2:2).  It  approximated  a  blackbody  absorber/radiator  when 
exposed  to  the  heated  samples.  On  the  other  hand,  the 
quartz  bell  jar  used  in  this  study  used  no  external  cooling. 
It  was  assumed  that  the  high  energy  radiation  from  the 
sample  was  either  transmitted  through  the  bell  jar  to  the 
surrounding  room  or  absorbed  by  the  jar  (2:9).  Other  high 
energy  sources  (i.e.  open  flames,  sunlight)  were  not  allowed 
to  be  present  during  experiment  runs.  Finally,  the  sample 
was  mounted  off-center  in  order  to  preclude  a  focusing 
effect  from  the  minimal  reflections  of  high-energy  radiation 
from  the  bell  jar  wall  (2:9). 

DC  Power  Supply 

Two  different  Intermagnetics  DC  power  supplies  were 
used.  Both  were  designed  for  high-current,  low-voltage 
electromagnet  applications.  Rated  at  75  amperes  maximum, 
the  first  supply  proved  insufficient  for .heating  the  Nb-1% 

Zr  sample  appreciably  due  to  the  relatively  low  resistivity 
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of  Nb-1%  Zr.  Therefore,  the  first  supply  was  replaced  with 
one  rated  at  180  A  maximum.  Each  contained  a  precision 
resistor  in  series  with  one  pole  of  the  output  lines.  Mea¬ 
suring  the  voltage  drop  across  the  resistor  permitted  calcu¬ 
lation  of  current  at  up  to  six  decimal  places  of  precision. 
Since  the  resistor  in  each  supply  was  calibrated  to  provide 
a  drop  of  100  mV  at  full  rated  current,  a  precision  digital 
voltmeter  was  used  as  a  handy  current  meter  by  displaying 
the  percentage  of  maximiam  current  (either  75  A  or  180  A) 
flowing  through  the  specimen. 

Due  to  the  relatively  high  current  involved,  power  was 
supplied  to  the  specimen  via  large  insulated  copper  cedsles 
of  00  gauge.  For  reasons  explained  in  Section  IV,  the 
experimental  procedure  required  a  reversal  of  current  during 
measurements  at  each  data  point.  Since  reversing  these 
large  cables  took  enough  time  to  allow  the  temperature  of 
the  sample  to  fall,  it  was  determined  that  mechanical 
switching  would  be  preferable.  Using  available  copper 
stock,  the  AFIT  Model  Shop  fabricated  a  double-pole,  double¬ 
throw  knife  switch  for  this  purpose.  The  knife  switch  was 
wired  into  the  circuit  so  that  it  reversed  the  polarity  when 
thrown.  In  order  to  ease  the  strain  of  sudden  open  circuit 
conditions  on  the  low- impedance  power  supplies  during  polar¬ 
ity  reversal,  a  copper  shorting  bar  was  used. 
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A  Hewlett  Packard  Model  3478A  Multimeter  with  six  digit 
display  was  used  for  taking  thermocouple  emf,  specimen 
voltage  drop,  and  precision  resistor  voltage  drop  measure¬ 
ments.  Twenty  times  more  sensitive  than  required  by  ASTM  C 
835,  the  meter  had  a  sensitivity  of  0.1  /iV.  Although  not 
exploited  due  to  time  constraints,  the  meter  possessed  a 
digital  data  output  capability  for  automatic  data  recording. 

Specimen  Holley 

As  previously  noted,  the  specimen  was  mounted  off- 
center  to  minimize  the  impact  of  normal  reflections  from  the 
bell  jar  wall.  Supporting  the  ends  of  the  specimen  were  two 
clamps  fashioned  from  woven  copper  braid,  flexible  enough  to 
allow  for  specimen  expansion  during  increasing  temperature. 
These  clamps  were  also  attached  to  standard  laboratory  cleunp 
stands  and  to  the  power  supply  cables.  As  shown  in  Figure 
4,  the  clamp  stands  also  served  to  support  the  thermocouple 
leads,  separating  them  from  the  specimen  and  from  each 
other. 

Thermocouples 

Also  referred  to  as  Type  s,  platinum/platinum- 10% 
rhodium  thermocouples  were  chosen  for  this  study  primarily 
for  the  temperature  range  they  span,  but  also  because  of 
their  ductility.  Suitable  positions  for  mounting  thermo¬ 
couples  were  scored  on  the  specimen  with  a  scalpel,  and  the 
distance  between  them  was  measured  under  a  microscope. 
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Figure  4.  Nb-1%  Zr  Specimen  Mounted  Inside  Bell  Jar 

Short  (few  cm)  lengths  of  3  mil  gauge  thermocouple  stock 
were  spot  welded  to  the  marked  positions  using  a  resistance 
welder  at  the  EG&G  Mound  Technologies,  Inc.  Facility.  A 
total  of  three  Pt,  Pt-10%  Rh  thermocouple  pairs  were  at¬ 
tached  at  approximately  4  mm  intervals.  Figure  5  shows  the 
relative  placement  of  the  thermocouples  and  power 
electrodes. 

Laser  welding  was  attempted  but  was  abandoned  due  to 
the  degree  of  difficulty  involved  in  working  with  niobium. 
Instead,  resistance  welds  using  copper  needles  proved  much 
more  effective.  In  order  to  preclude  the  occurrence  of 
oxidation  during  the  procedure,  a  procedure  developed  at 
Mound  was  employed  involving  the  direction  of  argon  through 
a  diffuser  stone  at  the  weld  at  a  flow  rate  of  6  cubic  feet 


per  hour  (9) . 


Figure  5.  Thermocouple  and  Electrode  Placement 

Although  Pt,  Pt-10%  Rh  thermocouples  have  response 
curves  published  up  to  1938  K,  a  eutectic  forms  at  the  weld 
between  Pt  and  Zr,  with  a  melting  point  of  1458  K.  Exceed¬ 
ing  1458  K  could  in  weld  failure  (11).  Other  eutectics 
exists  between  Nb  and  Pt,  Nb  and  Rh,  Zr  and  Rh,  etc.,  but 
the  lowest  melting  point  is  that  of  Pt  and  Zr.  For  that  - 
reason,  three  pairs  of  tungsten-2 6%  rhenium,  tungsten-5% 
rhenium  thermocouples  were  also  welded  to  the  specimen  in  an 
alternating  arrangement  in  anticipation  of  being  able  to 
continue  taking  data  at  temperatures  beyond  the  melting 
point  of  the  Pt-Zr  eutectic.  As  it  turned  out,  the  higher 
temperature  range  was  not  explored  due  to  power  supply 
constraints,  so  the  tungsten-rhenium  thermocouple  stiibs  were 
not  used. 


23 


At  this  point  it  is  necessary  to  mention  that  the 
apparatus  was  constructed  using  cold-worked  thermocouple 
material  which  was  not  annealed  prior  to  use.  Hindsight 
suggests  a  more  proper  technique  of  annealing  the  thermocou¬ 
ple  wire  stock  prior  to  building  the  feed  through  assemblies 
or  welding  the  pigtails  onto  the  sample.  In  addition,  the 
only  thermocouple  calibration  attempt  was  made  at  the  ice 
point.  Again,  a  more  prudent  procedure  would  probably  have 
included  some  method  of  obtaining  a  calibration  point  at 
high  temperature.  Both  these  problems  are  discussed  in  the 
results . 

A  thermocouple  feed-through  and  extension  assembly  was 
fabricated  by  the  AFIT  Model  Shop  using  a  plexiglass  flange, 
high-vacuum  epoxy,  and  10  mil  gauge  Pt,  Pt-10%  Rh  thermocou¬ 
ple  wire.  After  installing  the  assembly  into  a  port  on  the 
bell  jar  support  table,  the  10  mil  leads  were  threaded  into 
short  lengths  of  ceramic  insulation  and  were  resistance 
welded  to  the  3  mil  pigtails  protruding  from  the  specimen. 
Once  again,  argon  was  used  to  provide  an  inert  welding 
environment.  The  welds  between  the  two  gauges  of  wire  were 
then  cleaned  with  nitric  acid  to  remove  any  contamination 
acquired  during  assembly. 

An  aluminum-walled  insulated  ice  bath  was  fabricated  by 
the  AFIT  Model  Shop  for  use  as  a  thermocouple  temperature 
reference.  (The  ice  bath  was  necessary  due  to  the  fact  that 
the  published  thermocouple  response  tables  were  based  on  a 
reference  temperature  of  0*  C.)  Holes  were  drilled  into  the 
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lid  into  which  rubber- stoppered  glass  tubes  were  inserted. 
Each  tube  was  filled  with  a  small  amount  of  mercury  into 
which  a  single  thermocouple  lead,  connected  to  its  copper 
extension  wire,  was  inserted.  The  copper  extension  wires 
were  routed  through  a  multi-position  push-button  switch  to 
the  voltmeter.  For  ease  in  taking  measurements,  the  switch 
was  wired  so  that  each  of  the  three  thermocouples  could  be 
independently  selected.  Figure  6  depicts  the  wiring  ar¬ 
rangement  for  measuring  the  voltage  produced  by  thermocouple 
number  1 — the  switch  and  all  other  wires  are  omitted  from 
the  drawing  for  the  sake  of  clarity. 


Figure  6.  Wiring  Diagram  for  Thermocouple  Number  1. 
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In  addition,  two  other  positions  on  the  switch  allowed 
for  voltage  drop  measurements  across  the  sample  between  the 
top  (number  1)  ?.nd  middle  (number  2)  junctions  and  between 
the  middle  and  bottom  (number  3)  junctions  respectively. 
This  function  of  the  apparatus  is  demonstrated  in  Figure  7 
for  the  top  two  thermocouples.  Again,  only  the  active 
circuit  is  depicted. 


1 

2 

3 


Specimen  Ice  Bath 


Precision 

Voltmeter 


Figure  7.  Wiring  Diagram  for  Measuring  Voltage  Drop  Across 
Specimen  Between  Thermocouples  1  and  2 


Finally,  a  switch  position  was  reserved  for  measuring 
the  voltage  drop  across  the  power  supply's  precision  resis¬ 
tor,  shown  in  Figure  8.  Thus,  acquiring  data  at  each  point 
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was  relatively  simple — the  operator  had  only  to  record  the 
voltmeter  readout  corresponding  to  each  button  on  the 
switch. 


Figure  8 .  Wiring  Diagram  for  Measuring  Voltage  Drop  Across 

Precision  Resistor 
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lY.  Exp?riff9nt^l  PrQ<??durs 


Preparation 

The  apparatus  was  first  operated  using  a  thin,  flat 
specimen  of  oxidized  stainless  steel  (SS)  obtained  from  D. 
McElroy,  ORNL.  The  strip  was  approximately  20  cm  long,  1.25 
cm  vide  and  0.05  cm  thick.  It  ceune  prepared  with  three 
pairs  of  Pt,  Pt-10%  Rh  thermocouples  spot  welded  in  a  row 
near  the  center  of  one  side  about  2  cm  apart.  Configured 
essentially  the  same  as  the  Nb-1%  Zr  sample,  the  SS  specimen 
was  used  to  verify  the  experimental  procedure.  The  value 
for  the  emittance  of  this  specimen  had  been  previously 
measured  at  ORNL  using  an  apparatus  conforming  to  ASTM 
Standard  C  835.  The  dimensions  for  computing  surface  areas 
between  thermocouple  junctions  for  both  the  SS  and  the  Nb-1% 
Zr  specimens  had  been  previously  measured  under  a  microscope 
to  five  significant  digits  and  entered  into  the  notebook. 

Following  the  data  collection  runs  with  the  SS  speci¬ 
men,  the  apparatus  was  reconfigured  for  the  Nb-1%  Zr  tube. 
The  thermocouple  extensions  required  rewelding  to  the  Nb-1% 
Zr  specimen's  thermocouple  pigtails  and  the  electrodes  had 
to  be  reconfigured  to  accommodate  the  change  to  a  cylin¬ 
drical  geometry.  The  procedure  for  taking  data  was  essen¬ 
tially  the  same  for  both  specimens. 

Before  each  data  collection  run  was  made,  the  ice  bath 
was  prepared  using  a  mixture  of  de- ionized  water  and  cracked 
ice  made  from  the  same  water.  While  waiting  for  the  temper- 
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ature  of  the  bath  to  stabilize  at  the  freezing  point,  the 
following  pre-check  items  were  accomplished.  First,  a 
visual  inspection  was  performed  of  all  electrical  connec¬ 
tions  and  a  continuity  check  made  of  each  thermocouple  lead. 
Next,  the  vacuiom  gauge  reading  was  recorded.  For  safety 
concerns,  the  bell  jar  was  then  inspected  for  surface  cracks 
which  might  indicate  imminent  failure  of  the  quartz  under 
une^'en  heating  conditions.  The  power  supply  and  voltmeter 
were  turned  on  and  allowed  to  warm  up.  Finally,  when  all 
checks  were  complete,  the  power  supply  was  ramped  up  slowly 
to  allow  a  check  of  the  thermocouples  and  the  voltmeter.  At 
this  point,  an  arbitrary  determination  was  made  whether  to 
begin  taking  data  at  maximum  temperature  and  work  down  the 
scale  or  vice-versa.  Power  was  then  advanced  to  the  start¬ 
ing  point,  and  the  system  was  left  in  this  configuration 
until  the  temperature  stabilized.  Stabilization  was  deter¬ 
mined  to  be  reached  when  no  change  occurred  in  measurements 
of  thermocouple  emf  over  a  2  minute  interval. 

Data  Collection 

When  temperature  stabilization  was  reached,  the  system 
was  considered  to  exist  in  a  state  of  thermodynamic  equilib¬ 
rium.  A  check  of  the  ice  bath  was  then  made  to  ensure  the 
proper  reference  temperature  was  available  for  the  thermo¬ 
couples.  Data  collection  began  by  entering  the  following 
data  into  the  notebook; 
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1.  Time  of  Day 

2.  Polarity  (forward  or  reverse) 

3.  Current  as  percent  of  maximum,  read  directly  from 
voltmeter  (See  section  III.) 

4.  Emf  measured  from  thermocouple  1  in  V 

5.  Emf  measured  from  thermocouple  2  in  V 

6.  Emf  measured  from  thermocouple  3  in  V 

7.  Voltage  drop  between  thermocouples  1  and  2  in  mV 

8.  Voltage  drop  between  thermocouples  2  and  3  in  mV 

Following  collection  of  the  edaove  data,  the  power 

supply  was  short-circuited  and  the  knife  switch  thrown  to 
the  opposite  position  in  order  to  reverse  the  polarity  of 
the  current  applied  to  the  specimen.  As  soon  as  the  switch 
was  thrown,  the  shorting  bar  was  removed.  An  explanation 
for  the  polarity  reversal  and  the  short  circuit  procedure 
follows. 

The  purpose  for  reversing  polarity  was  to  provide 
another  emf  measurement  to  counteract  bias  error  introduced 
by  slightly  misaligned  thermocouple  wires.  The  mean  of  the 
two  opposite-polarity  measurements  was  taken  and  used  as  the 
"true"  emf  at  the  junction.  Each  reversal  was  accomplished 
after  first  shoirting  out  the  low-impedance  power  supply  to 
protect  it  from  frequent  full-load  to  open-circuit  transi¬ 
tions  and  vice  versa. 

When  polarity  reversal  was  accomplished,  the  sample  was 
allowed  to  return  to  equilibrium,  determined  by  noting  a 
steady  temperature.  Then,  the  data  collection  procedure  was 
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repeated,  except  for  the  polarity  reversal.  Thus  ended  the 
collection  at  the  initial  specimen  temperature. 

Next,  the  power  supply  was  adjusted  either  up  or  do%m, 
as  required,  by  a  predetermined  amount.  The  polarity  used 
in  the  last  half  of  the  prior  collection  run  was  used  in 
order  to  reduce  by  half  the  number  of  power  Interruptions  to 
the  sample.  After  allowing  the  system  to  reach  equilibrium, 
another  set  of  data  was  taken.  Then  the  polarity  was  re¬ 
versed,  and  the  process  repeated  until  the  predeteinnined 
ending  point  was  reached. 

At  various  intervals  during  the  procedure,  the  tempera¬ 
ture  of  the  ice  bath  was  checked,  and  ice  was  added  as 
necessary.  Bell  jar  pressure  was  continually  monitored. 

See  Appendix  B  for  sample  data.  Table  III  shows  sample 
data  using  the  oxidized  SS  specimen,  and  Table  IV  shows  the 
data  for  the  Nb-1%  Zr  specimen  before  visible  oxidation  set 
in. 


Data  Reduction 

National  Bureau  of  Standards  (NBS)  Circular  125  was 
consulted  to  obtain  thermocouple  response  curves  ( 17 ) .  The 
reference  listed  emf  potentials  for  Ft,  Pt-10%  Rh  thermo¬ 
couples  as  an  expanded  power  series  in  temperature  (See 
Appendix  A) .  An  iterative  solution  was  necessary  to 
calculate  temperatures  from  emf  data  measurements.  For  this 
effort,  TK  Solver  Plus,  Academic  Edition,  was  employed  (16) .~ 
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Once  the  temperatures  were  obtained,  the  mean  tempera¬ 
ture  was  calculated  for  each  pair  of  values  obtained  from 
thermocouples  1  and  2,  and  used  to  represent  the  temperature 
of  the  zone  between  the  two  thermocouples,  designated  the 
top  zone.  Likewise,  the  temperature  for  the  bottom  zone  was 
obtained  by  averaging  those  from  thermocouples  2  and  3. 

Next,  the  voltage  drop  between  thermocouples  1  and  2 
was  multiplied  by  the  current  flowing  through  the  sample  to 
obtain  the  amount  of  electric  power  being  dissipated  within 
the  top  zone.  Likewise,  the  power  for  the  bottom  zone  was 
calculated  using  thermocouples  2  and  3.  The  power,  tempera¬ 
ture,  and  surface  area  (previously  measured)  were  entered 
into  Eg  (13)  to  obtain  the  emittance.  This  was  done  for 
both  the  stainless  steel  and  the  Nb-1%  Zr  specimens.  During 
the  error  analysis,  it  was  determined  that  Tg  could  not  be 
ignored.  Tg  was  assumed  to  be  295  K  for  all  data  points. 
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Y.  Bsg^lta 


Introduction 

This  section  contains  an  analysis  of  the  results  of  the 
emittance  measurements  and  a  treatment  of  the  associated 
errors.  Results  for  the  stainless  steel  emittance  standard 
are  presented  first,  followed  by  the  results  for  the  fuel 
cladding  specimen.  The  method  used  for  estimating  the  error 
in  the  reported  emittance  values  is  taken  from  ASTM  C  835 
standard  and  is  sho%m  in  Appendix  C. 

Specimen  Emittance  ValVtgg 

Emittance  measurements  were  made  of  the  SS  strip  in 
this  study  from  ediout  780  K  up  to  about  945  K.  See  Appendix 
B,  Table  V  for  a  sample  of  initial  calculation  results. 
Figure  9  shows  an  initial  comparison  between  the  results  of 
this  method  and  those  from  R.  S.  Graves  and  D.  McElroy  of 
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Figure  9.  Initial  Results  Compared  to  ORNL,  Emittance  vs 

Temperature,  SS  Specimen 
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ORNL  (6) .  This  Study  yielded  an  initial  value  for  the 
emittance  of  eUsout  1.02  at  900  K  as  compared  to  a  linearly 
interpolated  value  of  0.88  for  the  ORNL  data.  A  linear 
interpolation  for  the  ORNL  data  was  determined  to  be  appro¬ 
priate  in  this  case  due  to  the  availability  of  only  two  data 
points. 

At  this  point,  it  was  apparent  that  the  thermocouple 
calibration  error  of  -34.8  K  at  room  temperature  was  preva¬ 
lent  at  higher  temperatures  and  could  not  be  ignored. 

Without  a  high  temperature  calibration  point,  it  was  assumed 
that  the  thermocouples  responded  similarly  throughout  the 
entire  temperature  range  used.  A  calibration  at  a  high 
temperature  was  not  performed  due  to  equipment  and  time 
constraints.  At  one  point,  an  optical  pyrometer  was  used  to 
examine  the  surface  of  the  heated  specimen.  Although  the 
pyrometer  was  useful  for  examining  the  relative  temperature 
along  the  surface  of  the  specimen,  no  actual  temperature 
calibration  was  possible  due  to  the  following  reasons. 

First,  the  emittance  of  the  object  being  obseiri^ed  was  an 
input  parameter  for  the  pyrometer  control  panel,  and  emit¬ 
tance  was  the  unknown  quantity  being  measured.  Second,  the 
absorption  characteristics  for  the  bell  jar  were  unknown  for 
the  spectral  band  at  which  the  pyrometer  was  most  sensitive. 

In  conclusion,  calibration  error  was  removed  from  the 
data  by  adding  34.8  K  to  all  temperatures.  Figure  10  shows 
one  data  set  from  Figure  9  after  applying  the  error  correc¬ 
tion.  The  set,  designated  "AFIT  22  Dec  88  (2)"  corresponds 
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Figure  10.  SS  Specimen  Amittance  Corrected  for  Thermocouple 

Calibration  Error 


to  Table  V  for  Figure  9  and  Table  VI  for  Figure  10.  Note 


that  Figure  10  also  includes  error  bars;  the  sources  for 


error  are  discussed  in  the  following  section. 


SS  Specimen  Error  Analysis 

According  to  T.  G.  Kollie  and  R.  S.  Graves,  cold  worked 
Type  S  thermocouples  can  introduce  a  maximum  error  of  -8  K 
(12) .  As  time  passed  during  experiment  runs,  the  ends  of 
the  thermocouples  nearest  the  specimen  annealed,  decreasing 
this  error  by  some  unknown  function  of  time,  temperature, 
and  proximity  to  the  specimen.  However,  this  error  was  part 
of  the  calibration  error  mentioned  earlier.  No  separate 
treatment  was  therefore  required. 

The  temperatures  obtained  from  the  three  thermocouples 
varied  by  about  4  K  maximum  at  any  given  data  point.  There¬ 
fore,  it  was  estimated  that  the  temperatures  reported  were 
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accurate  to  ±  2%.  The  uncertainty  in  the  emittance  for  the 
stainless  steel  specimen  typically  ranged  from  about  ±  3.6% 
to  ±  7.6%  (See  Appendix  B,  Table  VI  for  the  corrected  SS 
calculations  with  estimated  uncertainty.) 


Nb-1%  Zr  Specimen  Emittance  Values 

The  first  set  of  measurements  made  using  the  Nb-1%  Zr 
specimen  were  taken  to  represent  the  "virgin”  (unoxidized) 
condition.  The  values  ranged  from  0.159  ±  1.34%  at  913  K  to 
0.200  ±  1.25%  at  1091  K.  The  specimen  underwent  a  steady, 
visible  increase  in  surface  oxidation  from  one  run  to  the 
next,  increasing  the  emittance  in  the  process.  The  highest 
value  recorded  was  0.339  ±  1.63%  at  1269  K.  Figure  11 
demonstrates  the  upward  trend,  and  compares  the  results  with 
a  correlation  made  by  Senor  and  others  of  three  existing 
data  sets  (15:206). 
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TEMPERATURE  [K1 

Figure  11.  Emittance  vs  Temperature — Nb-1%  Zr 
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The  correlation  was  used  as  a  "predicted"  value  for  this 
study  and  is  labeled  as  such  in  Figure  11.  Appendix  A, 

Table  II  contains  the  formula  for  the  correlation — a  third 
degree  power  series  expansion  in  temperature.  It  should  be 
noted  that  the  surface  characteristics  of  the  specimens  from 
which  the  correlated  data  were  taken  were  known  to  be 
smooth.  Surface  oxidation  levels  were  not  addressed. 

See  Appendix  B,  Table  VII  for  sample  emittance  calcula¬ 
tions  and  estimated  uncertainty.  Table  VII  corresponds  to 
the  data  set  labeled  "10  Mar  89"  on  Figure  11. 

Nb-1%  Zr  Specimen  Error  Analysis 

Besides  the  two  sources  of  systematic  error  in  tempera¬ 
ture  measurement  noted  for  the  SS  specimen,  an  additional 
source  appeared  for  Nb-1%  Zr.  Being  an  alloy  of  refractory 
metal,  Nb-1%  Zr  has  a  very  high  tendency  to  oxidize  when 
heated,  even  at  a  vacuum  chamber  pressure  specified  by  ASTM 
C  835-82.  As  the  experiment  progressed,  the  surface  of  the 
Nb-1%  Zr  sample  steadily  increased  in  oxide  content,  driving 
up  the  emittance  as  discussed  earlier.  Due  to  the  increase 
in  emittance  as  an  unknown  function  of  time,  temperature, 
and  a  combination  of  the  two,  an  analysis  of  the  statistical 
error  or  repeatability  was  not  performed.  The  ASTM  standard 
quotes  a  three-sigma  bias  of  5%  for  the  method  in  general 
(2:10) . 

It  should  be  noted  that,  although  not  quantified,  there 
were  several  other  possible  sources  for  systematic  error. 
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First,  the  assumption  that  the  bell  jar  approximated  a  black 
body  edasorber  at  high  specimen  temperatures  by  transmitting 
the  high-energy  radiation  might  have  induced  additional 
error  from  reflections  off  the  protective  screen  surrounding 
the  jar.  In  addition,  each  thermocouple  veld  contact  area 
might  have  differed  enough  to  affect  emittance  calculations, 
despite  the  general  guidance  given  in  the  ASTM  standard. 


VI .  Summary  and  Recommendations 


Summary 

This  study  was  an  attempt  to  measure  the  total  hemi¬ 
spherical  emittance  of  a  sample  specimen  of  Nb-1%  Zr  fuel 
cladding  for  the  SP-100  space  nuclear  reactor  in  a  yacuum  at 
operating  temperatures.  The  purpose  was  to  inyestigate  the 
radiatiye  efficiency  of  the  cladding  to  aid  in  the  reactor 
coie  design.  Total  hemispherical  emittance  characterizes 
this  efficiency  and  affects  the  quantity  and  placement  of 
reactor  core  emergency  heat  sinks. 

Unoxidized  Nb-1%  Zr  emittance  yalues  ranged  from  0.159 
±  1.34%  at  913  K  to  0.200  ±  1.25%  at  1091  K.  The  highest 
yalue  obtained  for  the  oxidized  sample  was  0.339  ±  1.63%  at 
1269  K.  After  initial  calibration,  it  was  discoyered  that  a 
significant  portion  (about  25%)  of  the  calibration  error  was 
a  systematic  error  caused  by  improperly  prepared  thermo¬ 
couples.  Steady  oxidation  of  the  specimen  surface  affected 
the  results,  as  can  be  seen  in  Figure  11. 

Recommendations 

The  study  was  essentially  a  foundation  for  what  was 
perceiyed  as  a  continuing  effort.  Therefore,  seyeral  recom¬ 
mendations  are  offered  for  the  purpose  of  refining  the 
ayoidance  of  systematic  error  as  well  as  increasing  the 
scope  of  the  effort  to  improve  the  applicability  of  the 
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research  to  the  SP-100  program.  A  discussion  of  systematic 
error  avoidance  will  occur  first. 

To  begin  with,  any  samples  used  in  emittance  measure¬ 
ments  should  have  a  detailed  exeunination  and  appraisal  of 
the  sample's  surface,  to  include  grain  size  and  the  extent 
of  oxidation  or  other  surface  conteusinants .  Preferzdsly, 
this  should  be  done  by  a  metallurgist  or  someone  who  is 
knowledgeed>le  in  metallic  surface  features.  In  providing 
such  an  appraisal  with  the  emittance  experiment  results,  the 
scientific  community  can  compare  data  from  multiple  sources 
in  a  more  meaningful  manner.  In  addition,  a  detailed  sur¬ 
face  appraisal  would  permit  a  meaningful  comparison  of 
l2d>oratory  emittance  measurements  with  theoretical  expres¬ 
sions  of  emissivitv  as  a  function  of  resistivity  (1:294). 

For  validity,  such  equations  require  the  use  of  several 
assumptions  about  the  surface  conditions  of  the  specimen. 

In  addition,  further  work  on  this  alloy  should  occur 
only  in  a  high-vacuum,  clean-room  environment.  While  useful 
for  verifying  the  experimental  procedure,  the  vacuum  systems 
available  for  this  study  were  not  adequate  to  prevent  sur¬ 
face  contamination  of  the  sample.  Despite  the  fact  that  a 
Nb-1%  Zr  sample  is  likely  to  continue  to  outgas  at  high 
vacuum,  a  turbo-molecular  or  ion  pump  system  would  certainly 
be  preferable  to  an  oil  diffusion  pump,  for  instance. 

Another  source  of  error  noted  in  this  study  was  the  use 
of  cold-worked  Type  S  thermocouples.  Use  of  Type  S  stock 
should  include  an  annealing  period  prior  to  installation  in 
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the  apparatus.  In  addition  to  annealing,  the  thermocouples 
should  be  calibrated  at  the  high  end  of  the  range,  and  not 
just  at  the  ice  point.  The  melting  point  of  gold  (gold 
point)  is  one  candidate  for  such  a  calibration;  several 
threads  of  gold  could  be  welded  to  the  specimen  near  the 
zone  of  interest.  As  the  temperature  of  the  specimen  is 
increased,  the  gold  will  eventually  melt  and  fall  off.  At 
the  point  which  this  occurs,  the  thermocouple  potential 
should  be  noted,  and  the  calculated  temperature  compared  to 
the  melting  temperature  of  gold.  Of  course,  the  character¬ 
istics  of  any  gold-niobium  and  gold-zirconium  eutectics  must 
be  investigated  beforehand — this  technique  will  not  work  if 
the  eutectic  melting  point  lies  below  the  gold  point.  If 
the  process  could  be  repeated  with  other  high  melting  point 
materials  which  would  not  contaminate  the  bell  jar,  a  cali¬ 
bration  curve  could  be  constructed. 

Yet  another  source  of  possible  error  was  assuming  no 
loss  due  to  radiation  from  the  inner  surface  of  the  test 
zone.  Since  a  real  fuel  pin  would  be  closed  at  both  ends, 
future  research  should  be  done  with  closed  cylinders  to 
eliminate  any  possibility  of  radiation  escaping  through  the 
holes. 

Now  that  error  avoidance  has  been  covered,  the  remain¬ 
der  of  the  recommendations  will  be  focused  on  improving  the 
applicability  of  follow  on  research  to  the  SP-100  program. 
These  recommendations  include  increasing  the  temperature 
range  and  preparing  a  realistic  specimen  surface. 
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First  and  foremost,  the  upper  temperature  range  must  be 
explored.  While  data  at  1200  K  is  important,  accurate  data 
in  the  temperature  region  approaching  2000  K  is  more  criti¬ 
cal  in  terms  of  SP-100  safety  engineering  and  core  design. 
There  are  several  ways  to  achieve  the  higher  temperatures. 

One  way  of  increasing  the  temperature  of  the  specimen 
would  be  to  increase  its  resistance.  Since  the  sample 
length  is  limited  by  the  bell  jar  dimensions,  the  specimen 
could  be  split  longitudinally  and  a  smaller  cross  sectional 
area  would  then  present  a  higher  resistance  to  the  current 
flow.  Thus,  the  same  current  flow  would  produce  more  heat. 
Alternatively,  notches  could  be  cut  near  each  end,  serving 
the  additional  purpose  of  flattening  the  longitudinal  tem¬ 
perature  distribution  (2:10). 

Another  solution  is  the  use  of  a  larger  power  supply. 

In  this  case,  the  power  supply  cables  will  become  more  mas¬ 
sive  and  may  present  a  dilemma  when  designing  feed-throughs 
and  specimen  holders.  A  water-  or  nitrogen-cooled  jacket 
could  be  fitted  to  the  00  gauge  cables  used  in  this  study; 
care  should  be  taken  that  the  cooler  cables  do  not  produce  a 
temperature  gradient  in  the  area  of  interest  on  the  specimen 
surface  from  increased  conduction  heat  loss. 

A  recommended  approach  for  obtaining  a  more  realistic 
specimen  surface  would  include  heating  a  sample  of  cladding 
in  liquid  lithium.  If  the  thermocouple  pigtails  were  at¬ 
tached  first,  the  sample  could  be  moved  directly  from  the 
lithium  environment  into  the  sample  holder  inside  the  bell 
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jar.  Providing  that  the  lithium- treated  sample  would  not 
degrade  the  vacuum  to  an  unacceptable  level,  this  would 
allow  a  comparison  with  previous  work  in  order  to  character¬ 
ize  the  possible  effects  on  emittance  of  interactions  be¬ 
tween  the  coolant  and  the  fuel  cladding. 

Finally,  some  thoughts  are  offered  concerning  oxidation 
and  its  engineering  implications.  One  might  ask  why  a  fuel 
pin  (and  therefore  an  emittance  measurement  sample)  should 
not  be  pre-oxidized  in  order  to  obtain  the  highest  emittance 
possible.  The  reason  for  not  using  such  an  approach  at  this 
point  is  the  uncertainty  involved  in  what  the  core  environ¬ 
ment  would  do  to  the  oxidation  layer  over  time.  Consider  a 
possible  consequence  of  using  oxidized  samples  in  the  labo¬ 
ratory  for  obtaining  the  emittance  value  used  in  designing 
the  reactor  core.  If  the  actual  reactor  core's  environment 
eventually  removed  or  otherwise  adversely  affected  the 
cladding  oxidation  layer,  the  cladding  emittance  could 
decrease  below  the  minimum  level  required  for  removing  decay 
heat  properly. 
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APPENDIX  h 


Bangs 


Xsca 


-50  to 


630.74  *C 


630.74  to 


1064.43  'C 


5.3995782346  X  10® 
1.2519770000  X  10** 
-2.2448217997  X  10*® 
2.8452164949  X  10'® 
-2.2440584544  X  10*" 
8.5054166936  X  10*’® 


-2.9824481615  X  10* 
8.2375528221  x  10° 
1.6453909942  X  10*® 


Table  II.  Power  Series  Expansion  for  e  of  Nb-1%  Zr 

(15:206) 


400  to  2700  K 


6.39  X  10* 
4.98  X  10  ® 
3.62  X  10*® 
-7.28  X  10*’* 
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APPENDIX  fi:  Sample  Data  AZld  galgUlation  Rgg^ltg 


S2UBple  data  collected  for  both  the  oxidized  stainless 
steel  (SS)  emittance  standard  and  the  Nb-1%  Zr  specimen  are 
shown  as  Teddies  III  and  IV,  respectively.  TeUble  V  shows  an 
initial  look  at  raw  emittance  calculation  results  for  the  SS 
standard  prior  to  temperature  error  correction.  Table  VI 
shows  the  corrected  results  for  the  seune  specimen.  Table 
VII  shows  the  corrected  emittance  calculation  results  for 
the  unoxidized  Nb-1%  Zr  specimen. 
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emf-1 ,2.3R  emf  gen«ratod  during  'reverse*  polarity,  ±1E-04  V 

V-top-F,R  Voltage  drop  between  thermocouples  1-2,  *forward”  or  'reverse*  polarity,  ±1E--06  V 
V-bot-F,R  Voltage  drop  between  thermocouples  2  -  3,  ±1E-06V 
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Table  V.  Sar  ole  Initial  Calculation  Results- -Stainless  Steel  (Uncorrected) 
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Table  VI.  Sample  Calculation  Results — Stainless  Steel  (Corrected) 
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Appendix  C:  Estimation  of  Errors  and  Uncertainties 


Temperature  Calculation  Error 

As  discussed  in  Chapter  5,  the  thermocouples  exhibited 
a  known  calibration  error  of  ST  =  -34.8  K.  Approximately 
one- fourth  of  this  error,  or  -8  K,  was  due  to  improper 
preparation,  i.e.  using  cold-worked  instead  of  annealed 
thermocouple  stock.  Although  some  annealing  undoubtedly 
occurred  on  the  end  of  each  thermocouple  near  the  heated 
specimen,  the  maximum  error  of  -34.8  K  was  assumed. 

Tables  V  and  VI  compare  the  initial  results  with  the 
corrected  results  for  the  Stainless  steel  specimen.  Table  V 
shows  the  results  of  one  experimental  run  without  correcting 
for  the  temperature  error  mentioned  above.  Although  a 
temperature  dependence  was  demonstrated  mathematically,  the 
calculated  emittance  values  were  greater  than  unity,  and 
thus  physically  meaningless.  Table  VI  shows  the  same  data 
table  after  correcting  for  ST  by  adding  34.8  K  to  each 
calculated  value  of  T.  All  Nb-1%  Zr  data  is  shown 
corrected. 

Uncertainties 

To  estimate  the  uncertainty  in  the  calculated  emittance 
values,  the  method  outlined  in  ASTM  C  835  was  used  (2:10). 
First,  an  expression  was  developed  for  calculating  the 
uncertainty  associated  with  each  of  the  variables  in  Eg 
(13),  using  Gibbings'  text  "The  Systematic  Experiment” 
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(5:316).  Expressed  in  percent,  these  uncertainties  were 
calculated  at  each  data  point.  The  overall  uncertainty  in 
the  emittance  value  at  each  point  was  estimated  as 


+  o 


1/2 


(14) 


where  each  a  is  the  uncertainty  in  the  respective  term  in 
Eq  (13)  associated  with  its  subscript.  A  sample  calculation 
is  presented  below  for  the  first  data  point  in  Table  VII. 

First,  uncertainties  in  the  power  dissipated  term  shall 
be  discussed.  Since  a  wattmeter  was  not  used  in  this 
measurement,  the  power  must  be  expressed  in  terms  of  the 
values  measured.  This  is  given  by 


Pd  - 


(15) 


where 

Pp  =  Electric  power  dissipated  in  test  region 
Er  =  Voltage  drop  across  precision  resistor 
Rp  =  Resistance  of  precision  resistor 
Eg  =  Voltage  drop  across  test  region 

By  applying  the  chain  rule  of  differential  calculus  to  Eq 

(15),  one  can  observe  how  uncertainties  in  the  independent 

variables  will  affect  Pq.  Here  one  must  assume  that  the 

individual  tolerances  of  the  variables  involved  are 

equivalent  to  the  respective  standard  error  a.  If  OpD 

represents  the  standard  error  in  PD,  then 
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(16) 


An  uncertainty  for  Pj,  may  then  be  expressed  as  a  percent  of 

the  value  of  Pj,  at  each  data  point.  The  values  used  for 

estimating  the  uncertainty  in  P^  in  the  case  of  the  first 

data  point  in  Table  VII  are  as  follows: 

Pd  =  1.135  W 

Es  =  15.291  mV 

Er  =  39.904  mV 

Rp  =  5.55E-04  ohms 

Ce,  =  lE-06  V 

Ubp  =  lE-06  V 

=  1% 

Next,  the  uncertainty  in  the  surface  area  measurement 
shall  be  discussed.  The  formula  for  the  surface  area  of  the 
test  region  is  A  =  jrdL  ,  where  d  is  the  outside  diameter 
of  the  tubing  and  L  is  the  length  of  the  test  region.  Given 
that  the  same  instrument  was  used  to  measure  the  lengths  of 
the  test  regions  and  the  diameter  of  the  tubing,  the 
uncertainty  in  area  was  estimated  as  twice  the  measurement 
uncertainty  in  one  dimension.  (This  estimation  held  true 
for  the  stainless  steel  specimen  as  well,  since  the  same 
instrument  was  used  to  measure  the  length  and  width  of  the 
test  regions.)  Since  the  uncertainty  in  one  dimension  was 
estimated  as  1.0  X  10'*  inches  (2.5  X  10"’  m)  ,  the  estimated 
maximum  uncertainty  in  surface  area  was  5.0  X  10"’  m*.  Using 
the  smaller  of  the  two  test  regions  for  the  Nb-1%  Zr 
specimen  (area  of  1.7788  X  10‘*  m^)  ,  the  estimated  maximum 
uncertainty  in  surface  area  was  given  by  /  A  =  0.28  % 
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Finally,  the  uncertainty  in  the  temperature  term  was 
estimated  by  considering  both  temperature  averaging  and 
voltmeter  tolerance.  First,  the  temperatures  reported  were 
in  fact  average  test  region  temperatures.  Although  the  test 
region  was  supposed  to  exhibit  a  constant  surface 
temperature,  there  was  a  variable  difference  noted  between 
adjacent  thermocouples.  An  examination  of  the  data  revealed 
a  maximum  difference  of  about  4  K  at  any  given  data  point. 
Second,  the  voltmeter  tolerance  for  the  range  used  in 
measuring  thermocouple  emf  was  1  X  10"*  V,  resulting  in  an 
insignificant  (far  less  than  1  K)  effect  on  temperature. 
Therefore,  with  a  maximum  4  K  spread  across  a  test  region, 
it  was  determined  that  a  region's  temperature  could  not  be 
measured  any  more  accurately  than  ±  2  K.  This  uncertainty 
can  be  expressed  as  a  percent  at  each  data  point.  However, 
the  ASTM  standard  uses  a  format  of  reporting  the  uncertainty 
in  the  [T*  -  T*o]  term  of  Eg  (14),  which  is  four  times  the 
uncertainty  in  T. 

The  results  of  a  sample  calculation  for  estimated 
uncertainty  is  shown  below  in  Table  VIII. 

Table  VIII.  Sample  Uncertainty  Estimation 


Emittance 

0.159 

Temperature ,  K 

913 

a,,  ±% 

0.97 

‘^[T^-To*]  ' 

0.88 

0.28 

«^PD/ 

0.98 

(First  data  point 

in  Table  VII) 
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